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Abstract 



An efficient method to calculate tadpole diagrams is proposed. Its capability is 
demonstrated by analytically evaluating two four-loop tadpole diagrams of current 
interest in the literature, including their O(e) terms in D = 4 — 2e space-time 
^ ■ dimensions. 
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<D . 1 Introduction 

> 

^ . Calculations of higher-order corrections are very important for precision tests 

of the Standard Model in present and future high-energy-physics experiments. 
The complexity of such calculations and the final results strongly increases 
with the number of quantum loops considered, and one rapidly reaches the 
limits of the present realm of possibility when one attempts to exactly account 
for all mass scales of a given problem, already at the two- loop level. To simplify 
the calculations and also the final expressions, various types of expansions 
were proposed during the last couple of years (see Ref. [1,2,3] and references 
cited therein). These approaches usually provide a possibility to reduce the 
problem of evaluating complicated Feynman integrals to the calculation of 
tadpoles T m (ai, a 2 ) and loops with massless propagators L q (ai, a 2 ), defined 
in Eqs. (2) and (4), respectively, which have simple representations. 
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As a first step, let us to introduce some definitions to be used below. All the cal- 
culations are performed in Euclidean momentum space of dimension D = 4— 2e 
using dimensional regularisation with 't Hooft mass scale //. For the Euclidean 
integral measure, we use the short-hand notation Dk = n' i ~ D d D k/n D l 2 . The 
dotted and solid lines of any diagram correspond to massless and massive 
Euclidean propagators, represented graphically as 

_*_ = = , a , (1) 

(k 2 ) a ' (k 2 + m 2 ) a m 2 y } 

where a and m denote the index and mass of the considered line, respectively. 
Unless stated otherwise, all solid lines have the same mass m. Lines with index 
1 and mass m are not marked. 

Let us now discuss the rules used in our calculation. Firstly, massive tadpoles 
To tm (eti, 0.2) are integrated according to the identity 



R(ai,a 2 ) 



where 



(2) 



%)0]| = I(M^M (3) 
V ' ; T(a 2 )V(D/2) v ; 



and r is Euler's Gamma function. Secondly, massless loops L 9 (q; 1 ,q; 2 ) with 
the external momentum q are integrated according to the identity 



L q (a 1 ,a 2 ) = J 



Dk 



(k 2 )^(k-q) a2 " q 




«2 



A( ai ,a 2 ) , (4) 



where 



o(ai)o(a 2 ) , , V{D/2-a) 

A(a 1 ,a 2 ) = — ■ — — , a a = — — . 5 

a(ai + a 2 — D/2) r(a) 



Recently, two examples of four-loop tadpoles were calculated in Ref. [4]. Cer- 
tain 0(e°) parts of these results, denoted as N 10 and A^o, were presented there 
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only numerically. 2 Also the 0(e) terms beyond Niq and N 2 o are of current 
interest in the literature [6] and so far only known numerically. In conjunction 
with the results of Ref. [7], they allow us to derive analytic results for several 
four-loop master integrals which are indispensable for modern calculations at 
this level of accuracy [8]. 

The subject of this letter is to advertise a powerful technique to evaluate 
such tadpoles analytically. The simplifications encountered in the examples at 
hand suggest that this technique might also be useful for more complicated 
diagrams. Therefore, we wish to introduce it to the interested reader. 

The content of this letter is as follows. Section 2 explains the calculational 
technique in general. Sections 3 and 4 describe its application to the two 
integrals iV 10 and N 20 of Ref. [4], respectively, also including their 0(e) terms. 
A summary is given in Section 5. 



2 Technique 



The core of the technique is a master formula to represent a loop with two 
massive propagators as an integral whose integrand contains a new propagator 
with a mass that depends on the variable of integration. This formula has the 
following form: 



Oil, ™1 



CK2 5 W12 



_ r(c*i + a 2 - D/2) 

r(ai)r(a 2 ) 

j ds ai +a 2 -D/2 

X J (1 - s )«i+l-D/2 s a 2 +l-D/2 q * m 2 m 2 — • W 
T=7+ -t 

It was introduced in Euclidean and Minkowski spaces in Refs. [9] and [10], 
respectively. Here, we are working in Euclidean space and will, thus, fol- 
low Ref. [9]. Some recent applications in Minkowski space may be found 
in Ref. [11]. Special cases of Eq. (6) were considered in Ref. [12], where a 
differential-equation method [12,13] was introduced. The latter only requires 
that quite simple diagrams are calculated directly. The results for more com- 
plicated diagrams may be reconstructed by integrating the results for simpler 
diagrams over external parameters. In many cases, the results for complicated 
diagrams may be obtained by integrating the one- loop tadpoles of Eq. (2). 



2 During the preparation of this article, we were informed by K.G. Chetyrkin about 
a paper [5] which also contains analytic results for Niq and N 2 o- Our results are in 
full agreement with those of Ref. [5] . 
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The method is now very popular for the calculation of complicated Feynman 
integrals; for recent articles, see Ref. [14] and references cited therein. 

Here, we shall follow a similar strategy. Applying Eq. (6), we shall represent 
the results of Ref. [4] as integrals over one-loop tadpoles, which in turn contain 
propagators with masses that depend on the variable of integration. The case 
of N w will be considered in detail. 

We would like to note that Eq. (6) can be applied successfully to diagrams con- 
taining one-loop self-energy sub-diagrams. Such sub-diagrams are frequently 
generated by the application of the integration-by-parts technique [15]. Such 
an application [12,13] leads to a differential equation for the original diagram 
(in exceptional cases, the equation degenerates to an algebraic relation) with 
an inhomogeneous term that depends on less complicated diagrams usually 
containing one-loop self-energy sub-diagrams (see, for example, Refs. [9,10]). 



3 First example 



As the first example, we consider the case of N\q in Ref. [4], which may be 
represented graphically as 3 



m 





T 4 (l + e) 



(m 



2\4e 



(7) 



where N w (e) = [iV 10 +M 10 e+O(e 2 )]/(l-e). For e = 0, we recover iV 10 = iV 10 (0). 
Applying Eq. (4) to the massless loop in the left diagram and Eq. (6) with 
q;i = 1 and a2 = 2 to the interior massive loops and to the loops with exterior 
propagators, we easily obtain the following representation for iVio(e): 



N 10 (e) 



m 
e 



1 



m 



2c 



■Mo) 



r 2 (i - e] 

r(l-2e 



-■Mi) 



where 

j . (1 — e)(m 2 ) 4e } ds\ f ds 2 m , , s 



3 The graphic representations of N\q and N20 are adopted from Ref. [4] . A propa- 
gator with a point is equal to one with the index 2. 
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Mf = m 2 /r]i, rji = Sj(l - s^, and T Mi ,m 2 ,m 3 («i, a 2 , a 3 ), where the index a { 
belongs to the mass M iy is the one-loop tadpole involving three massive prop- 
agators, 



^Mi,M2,M 3 («l, ^3) — 



Cx.2 /\ a 3 




(10) 



Expanding the propagator {k 2 + M|) ° 2 in To i Mi,m 2 (q ; 3, &i, « 2 ) as 



E 



T(n + a 2 ) /i™2 



(P + M|) Q 2 ^ n!r(a 2 ) (A; 2 + Mf) n+1+e ' 



11^ 



where ^ = Mf — M 2 , and using Eq. (2) for the tadpole T 0i Mi(aii, ^2), we 



obtain 



rr , x R{a 3 ,a 1 + a 2 ) ,_ , , , , 

Jo,Mi,M 2 (a3, "l, «2j = ^ M2 ^ D/2 2^i{a - D/2, a 2 ; a x + a 2 ; x), (12) 



where a = a\ + a 2 + 0:3, x = 1 — rr, rr = M|/M 2 = i]i/i] 2 , and 2 i*\ denotes a 
hyper-geometric function [16]. Setting ctj = 1 + a^e and expanding the hyper- 
geometric function, we have 



7o,Mi,m 2 (1 + a 3 e, 1 + aie, 1 + a 2 e) 



T(l + e) 



(M 2 ) 1+ ( 5+1 ) e x(l 



lnx 



+ e 



a 2 + a 3 + 1 2 . j ■ /_\ 
In 1- (ai + a 2 ) Li 2 {x) 



(a 2 + a 3 + l) 2 3 



6 



ln d x 



- ((a 3 + l)(a + 1) - l)C(2) lnx + a 1 (a 3 + l)lnxLi 2 ( 
+ (a 1 + a 2 ) 2 Li 3 (x)-(a 2 (a + l)-a 1 (a 3 + l))Si,2(3 T ) j + 0(e 3 ), (13) 



where a = ai + a 2 + a 3 and 



/_1 \n+m-l „ J 

Sn, m (^)= J n , . / ^ln^ln^l-xy), 
(n — l)!m! J y 

Li n (x) =Sn-l,l( a; ) 



(14) 



denote the generalised and ordinary poly-logarithms (see, for example, Ref. [17]), 
respectively. Thus, we obtain 
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N- 



(e) = / — - ( — — ^ { — In 2 x + In x In r/i 

J Si7?f J s 2 ?7? a; I 2 



10 



+ e 



19 



ln 3 x + 4C(2)lna;-lnxLi 2 (x) - (3 In 2 x - 2 Li 2 (x)) In 771 



} + 0(, 2 ). 

Exploiting the Si <-> s 2 symmetry, we find 



+ - In x In 771 



(15) 



1 1 



iV- 



( e ) = ^ / ds ± / — (lnxln(r?ir? 2 ) + e 
8 J J 779 — 771 I 



10 







^2 - ?7i 



61 
12 



In x 



+ 8C(2) lux + (Li 2 (3f) - Li 2 ln(7 7l 7 72 ) 



+ - lnxln 2 (77i?7 2 ) 



} + o( e 2 



(16) 



We first concentrate on the O(e ) term AT 10 . Using the standard replacement 
Sj = (1 + £i)/2, we obtain 



S1 s/ 

if. f 



'i n 2 LJi _ in 2 LJI 



(17) 





Exploiting the £1 <-> £ 2 symmetry, we find 



iV 10 = / dfc [ln 2 (l - £ 2 ) - 21n41n(l - £ 2 )] | ^ 



e 2 



= 4 / f 1 111 rrf I ln2(1 " e2) - 2 1114 111(1 " 



(18) 



Evaluating the r.h.s. of Eq. (18), we obtain 



N- 



17tt 4 7 



10 : 



+ -C(3)ln2-4Si, 3 (-l) 



720 2 
49tt 4 1 

"720" ~~ 2 &4 ' 



(19) 
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where 



6 4 = --(tt 2 - In 2 2) In 2 2 + 8 Li 4 f- ) . (20) 



Equation (19) coincides with the result obtained in Ref. [5]. 

The evaluation of the 0(e) term M\o is rather tedious and cannot be described 
in this brief communication. Here, we merely list the result, which reads 

14Q7T 4 279 

M 10 = -^ r ln2 + ^C(5)-76 5 , (21) 



where 



b 5 = — (5tt 2 - 3 In 2 2) In 3 2 + 8 Li 5 (- ) . (22) 



We note in passing that also the other master integrals presented in Ref. [5] can 
be written in a more compact form if the combinations 64 and 65 of Eqs. (20) 
and (22), respectively, are introduced. 



4 Second example 




Let us consider the second example of Ref. [4], which may be graphically 
represented as 



m2 l X * = w* (f '' (23) 



where A^ 20 (e) = [A^ 20 +M 20 e+O(e 2 )]/(l-e). For e = 0, we recover N 20 = iV 20 (0). 
Applying Eq. (6) with a\ — 1 and a 2 = 2 to the two interior massive loops 
and to the loops involving exterior propagators, we easily obtain the following 
representation for iV 2 o(e): 

(1 -e)(m 2 ) 4e f dsi f ds 2 f ds 3 
N 20 {e) = — / / / -T Mi ,m 2 mA 1 + e, 1 + e, 1 + e).(24) 



T(l + e) J sirjf J S2V2 { s 3 ^l 



Expanding the propagators (k 2 + Mf) ° 2 and (k 2 + Mf) 03 in 
^Mi,M 2 ,Af 3 (Q!i) «2, CK3) as in Eq. (11) and using Eq. (2) for the tadpole To,Mi(0, a 2 ), 



7 



we obtain 



( ,™r{n + a 2 )™r(l + a 3 )r{l + n + a-D/2) 

^M!,M2,m 3 «i, "2, a 3 J - 2^ — TFT — 7~ "7TF7 — \ F77~i I - =\ 

^ n!r(a 2 ) Jto '' r ( a 3j r(/ + n + a) 



^12^12, 



(M?) n+,+a - £ V 2 
#(0,a 



( M 2)a-D/2 



Fi(a - D/2,a2,a 3 ;a;x,fj, 1 3/M 1 ), 



where i*\ denotes an Appel hyper-geometric function [16]. Expanding the 
hyper-geometric function, we have 



T 



M 1 ,M 2 ,M ; 



3 (1 + aie, 1 + a 2 e, 1 + a 3 e) 



1 



r(i 



X 



M 2 M 2 



'2 

A*12 



a 2 + 1 M? 

2 /i!2 



(M 2 )( s+1 ) e /i 32 (l-e) 
M| /_M 2 



m —pr + a ai - a 3 Li 2 



2 . 

+ a 3 Li 2 

+ (2<->3). 



Li 2 



^21 
M 2 



Lio 



/^23 

M 2 



+ 0(6 2 



A*12 

M 2 



(25) 



As in the previous section, we describe the evaluation of the O(e ) term in 
some detail, while we merely list the final result for the 0(e) term for reasons 
of space. At 0(e°), the right-hand side of the Eq. (24) contains the tadpole 
T MltM2tM3 (l, 1, 1), which can be calculated in an essentially simpler way. In- 
deed, using 

1 1/1 1 

(P + M 2 )(A; 2 + M 2 ) ~ \k 2 + M? ~ k 2 + M 2 




and Eq. (2), we find 



7mi,m 2 (1, 1) 



r(c-i) 



A*2i 



(M 2 ) 1 - £ -(M 



1 r(l + e) / 1 M 2 M 2 \ 

— ^ M- + ^ln^- + (M 1 ^M 2 ) 

fi 2e 1 - e \e 1121 At 2 / 

+ 0(e), 

pMi,Af 2 (l> 1) — 7mi,M 3 (1) 1)] 

1^32 

M? M? 

In -4 + (Mi <-> M 2 ) + (Mi <-> M 3 ) 

A*i2A*i3 At 



(27) 
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+ 0(e). 



(28) 



Notice that the right-hand side of Eq. (28) is m 2 and fi 2 independent. More- 
over, it coincides with Eqs. (24) and (25) at O(e ). Every term in Eqs. (27) 
and (28) contains only one logarithm depending on the variable rji. Thus, the 
integrals Jq dsi can be evaluated at the end of the calculation (see Eq. (18), 
for a similar procedure). After some algebra, we obtain 

N20 = ^C(3), (29) 



which agrees with Ref. [5]. 

Lack of space prevents us from going into details with our derivation of the 
0(e) term M 20 , and we merely list our result, 

47T 4 

M 20 = -— - 27C(3) + 36 4 , (30) 
10 



where 64 is defined in Eq. (20). 



5 Conclusion 



We demonstrated the usefulness of a simple technique for the analytic evalua- 
tion of four-loop tadpoles, which allows one to obtain the results by integrating 
one-loop tadpoles with masses that depend on integration variables. By means 
of this method, we calculated two Feynman diagrams which were presented in 
Ref. [4] in numerical form. Our results are in full agreement with the analytic 
results obtained just recently in Ref. [5]. 

We note that the integral Niq was found in Ref. [5] with the aid of the PSLQ 
program [18], which is able to reconstruct the rational-number coefficients mul- 
tiplying a given set of transcendental numbers from a high-precision numerical 
result, i.e. the algebraic structure of the result must be known a priori. Such 
structures are often more complicated, and they are in general not known from 
other considerations (see, for example, Ref. [19]), so that the PSLQ program 
cannot be applied. The method presented here does not suffer from such a 
limitation and should be applicable also for more complicated tadpoles de- 
pending on one [5] or two [20] nonzero masses. Moreover, it can be applied 
in combination with the PSLQ program. In fact, by evaluating some part of 
a complicated Feynman diagram by means of this method, the underlying 
set of transcendental numbers may emerge and can then be injected into the 
reconstruction of the full analytic result with the help of the PSLQ program. 
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